Isotropic-Nematic Transition in Liquid-Crystalline Elastomers: 
Lattice Model with Quenched Disorder 
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When liquid-crystalline elastomers pass through the isotropic-nematic transition, the orientational 
order parameter and the elastic strain vary rapidly but smoothly, without the expected first-order 
discontinuity. This broadening of the phase transition is an important issue for applications of 
liquid-crystalline elastomers as actuators or artificial muscles. To understand this behavior, we 
develop a lattice model of liquid-crystalline elastomers, with local directors coupled to a global 
strain variable. In this model, we can consider either random-bond disorder (representing chemical 
heterogeneity) or random- field disorder (representing heterogeneous local stresses). Monte Carlo 
simulations show that both types of disorder cause the first-order isotropic-nematic transition to 
broaden into a smooth crossover, consistent with the experiments. For random-field disorder, the 
smooth crossover into an ordered state can be attributed to the long-range elastic interaction. 



I. INTRODUCTION 

Liquid-crystalline elastomers are complex materials 
consisting of crosslinked polymer networks covalently 
bonded to long, rigid, liquid-crystalline units 0,0,0] • Be- 
cause of this unusual structure, they combine the elastic 
properties of rubbers with the anisotropy of liquid crys- 
tals. Any distortion of the polymer network affects the 
orientational order of the liquid crystal, and, likewise, any 
change in the magnitude or direction of the orientational 
order influences the shape of the elastomer. These mate- 
rials have a low-temperature nematic phase, with long- 
range orientational order in the liquid-crystalline units, 
and a high-temperature isotropic phase, with no long- 
range orientational order. Near the isotropic-nematic 
transition, a small change in temperature induces a large 
change in the orientational order, which causes the elas- 
tomer to extend or contract substantially. This ther- 
mally induced extension and contraction enables liquid- 
crystalline elastomers to be used as actuators or artificial 
muscles HIE El 

A key issue for basic and applied research on liquid- 
crystalline elastomers is understanding the isotropic- 
nematic transition. In conventional liquid crystals of 
small molecules, this is a first-order transition, with a 
discontinuity in the magnitude of the orientational or- 
der as a function of temperature. By contrast, exper- 
iments on liquid-crystalline elastomers show that both 
the orientational order and parameter and the elastic 
strain vary rapidly but smoothly across this transition, 
with no first-order discontinuity 0, IE IE 0, IE H EJ- 
Apparently the experimental behavior is neither a first- 
nor a second-order transition, but rather a nonsingular 
crossover between the isotropic and nematic phases. Al- 
though the transition is sharp enough for applications, it 
is puzzling from the theoretical point of view. We would 
like to explain the broadening of this phase transition in 
liquid-crystalline elastomers, compared with the analo- 
gous transition in conventional liquid crystals. 

In a previous paper |Tl|. we considered two possible 



explanations for this broadening. The first possibility 
is that some aligning stress shifts the transition past a 
mechanical critical point [h2| . like a liquid-gas transi- 
tion at high pressure. This aligning stress might arise 
from an applied tensile stress on the sample, or from 
an anisotropic internal stress due to crosslinking an elas- 
tomer in the nematic phase. The second possibility is 
that the transition is broadened by heterogeneity in the 
elastomer. In particular, we considered heterogeneity in 
the local isotropic-nematic transition temperature. 

To assess these two possible explanations, we mea- 
sured the elastomer strain as a function of temperature 
over a range of applied tensile stress. We used elas- 
tomer samples crosslinked in the nematic phase, which 
should have a large anisotropic internal stress imprinted 
by the crosslinking process, and samples crosslinked in 
the isotropic phase, which should not have an anisotropic 
internal stress. By analyzing the experimental data, we 
found three indications that the broadening of the phase 
transition is caused by heterogeneity. First, the slope of 
strain vs. temperature at the transition does not depend 
sensitively on applied stress, in contrast with the predic- 
tion for homogeneous elastomers. Second, the broaden- 
ing occurs even for samples crosslinked in the isotropic 
phase. Third, the data for strain vs. temperature could 
not be fit well by the predictions of Landau theory for 
homogeneous elastomers, but could be fit much better by 
the homogeneous theory convolved with a heterogeneous 
distribution of transition temperatures. 

Although our previous study showed the importance 
of heterogeneity for the isotropic-nematic transition in 
liquid-crystalline elastomers, this study still leaves two 
open questions. The first issue is the distribution of local 
strains. The previous theory considered an average over 
local regions with different transition temperatures. At 
any given temperature, these regions have different local 
nematic order parameters and different local strains. It 
is not clear how regions with different local strains can fit 
together. The second issue is the type of heterogeneity. 
The previous theory considered a distribution of the lo- 
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FIG. 1: Schematic illustration of the lattice model for the 
isotropic-nematic transition in liquid-crystalline elastomers. 
In the isotropic phase, the directors are disordered, and there 
is no strain. In the nematic phase, the directors are ordered 
along one axis, and the material is extended with strain along 
that axis. 



cal transition temperature, which could arise from chem- 
ical heterogeneity in an elastomer. This type of hetero- 
geneity would be regarded theoretically as random-bond 
disorder. It is not the only possible type of heterogene- 
ity. Another possibility is a distribution of local stresses, 
which could arise from local orientational order in differ- 
ent directions at the time of crosslinking. This possibility 
would be regarded theoretically as random-field disorder. 
Several recent papers have considered random-field disor- 
der in liquid-crystalline elastomers [TiilT^llHll^ll^ll8j . 
These studies have shown that random fields strongly 
affect mechanical properties and correlation functions 
in the low-temperature nematic phase. However, they 
have not made predictions for the effects of random-bond 
or random-field disorder on the temperature-dependent 
isotropic-nematic transition. 

The purpose of this paper is to develop a lattice model 
for liquid-crystalline elastomers, which addresses these 
theoretical questions about the isotropic-nematic tran- 
sition. In this model, we explicitly consider both ori- 
entational order and elastic strain, as shown in Fig. 1. 
Orientational order is described by a nematic director, 
which is defined on each site of a three-dimensional (3D) 
lattice, as in the Lebwohl-Lasher model of liquid crys- 
tals . By contrast, elastic strain is defined by a single 
global lattice distortion variable. Because the strain is 
assumed to be uniform, there is no problem of fitting to- 
gether regions with different strains. In this model, we 
can consider either random-bond or random-field disor- 
der. Random-bond disorder enters the model as a vari- 
ation in the strength of the local coupling constant be- 
tween local directors on neighboring lattice sites, which 
implies a variation in the local transition temperature. 
By comparison, random-field disorder enters the model 
as a variation in the direction of an aligning field that 
acts on each local director, which implies a local stress 
on the elastomer. Both types of disorder are quenched, 
meaning that they are fixed and cannot evolve toward 
equilibrium. 

Our study leads to four main results. First, we derive 
a new theoretical formalism for liquid-crystalline elas- 



tomers, which translates the Warner- Terentjev neoclas- 
sical rubber elasticity 0- [H Q into a lattice Hamilto- 
nian for interacting directors and strain. Second, we use 
Monte Carlo simulations to determine the orientational 
order parameter and elastic strain of a homogeneous elas- 
tomer as a function of temperature and applied stress. 
These simulations show that the homogeneous elastomer 
has a mechanical critical point, as expected from Landau 
theory. Third, we simulate an elastomer with random- 
bond disorder, and find a broadening of the isotropic- 
nematic transition in both the orientational order and 
the strain. This result is consistent with experimental 
data for liquid-crystalline elastomers. Fourth, we sim- 
ulate an elastomer with random-field disorder, and find 
that the transition is also broadened in this case, con- 
sistent with experiments, provided that the random field 
strength is in the right range. The result for the random- 
field system is surprising, because random fields generally 
destroy long-range order rather than broadening a tran- 
sition to an ordered state. We attribute this result to an 
effective long-range interaction mediated by the elastic 
strain. 

The plan of this paper is as follows. In Sec. II we 
work out the theoretical formalism, leading to an explicit 
lattice Hamiltonian that can be simulated. In Sec. Ill 
we present the simulations for homogeneous, random- 
bond, and random-field elastomers, and give numerical 
results for the isotropic-nematic transition in each case. 
In Sec. IV we discuss these numerical results, and com- 
pare them with experiments and with other theoretical 
studies of quenched disorder. 



II. MODEL 

In order to simulate the isotropic-nematic transition 
in liquid-crystalline elastomers, we need a mathematical 
model for the interacting orientational and elastic degrees 
of freedom. In this model, we define a local nematic di- 
rector rti on each site i of a 3D cubic lattice. As in a 
conventional liquid crystal, the director rii is equivalent 
to — rij. The directors interact with a global lattice dis- 
tortion tensor A, which represents the overall shape of 
the sample. A schematic view of the directors and lat- 
tice distortion is shown in Fig. 1. 

The Hamiltonian for this lattice model can be written 
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^ ^ -^interaction ( 



n h flj) + ^ -^elastic (rij, A). (1) 



The first term in Eq. (JTJ is an interaction that favors 
alignment of the directors on neighboring lattice sites i 
and j. As in the Lebwohl-Lasher model, this interaction 
can be written explicitly as 



-^interaction ij^i t ) Jiji^i ' T^j 



(2) 



where Jy > is the local bond strength, which may be 
either uniform or disordered. The second term in Eq. 
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is an elastic term that couples the director orientation at 
lattice site i with the shape of the polymer chains, which 
are determined by the overall lattice distortion tensor A. 
This term should favor alignment of the directors along 
an orientation determined by the lattice distortion, and, 
conversely, favor a lattice distortion in an orientation de- 
termined by the directors. 

To develop an explicit expression for this elastic term, 
we use an argument based on the neoclassical rubber elas- 
ticity of Warner and Terentjev 0,0,0. In their theory, 
they derive the general formula 



p _ A 

r elastic ^ 



Tr (in ■ A T ■ £~ Y ■ A 



In 



det£ 



(3) 



where fi is the shear modulus, A is the lattice distortion 
tensor, £ is the shape tensor of the polymer chains, £o 
is the shape tensor at the time of crosslinking, and a is 
the average polymer step length. They further work out 
a specific expression for the shape tensor in the freely 
jointed chain model. That specific expression is not ap- 
propriate for a lattice Hamiltonian, because it is based 
on a global average over a system with imperfect nematic 
order. By contrast, in a lattice model, there is a director 
on each lattice site, which represents perfect local orienta- 
tional order along rii at site i. Any imperfect long-range 
order must emerge spontaneously from simulations of the 
interacting system. Hence, we should construct a local 
shape tensor £j for site i, derived from the director rij, 
which will enter into the trace formula of Eq. © . 

To construct the local shape tensor £i, we first con- 
sider a coordinate system aligned along the director n,, 
in which £i is diagonal. In that coordinate system, we 
have 



(4) 



where £± and £\\ are the anisotropic polymer step lengths 
favored by the local mesogenic unit. In a general coordi- 
nate system, the shape tensor components become 




(5) 



A similar argument gives an explicit expression for the 
lattice distortion tensor A. Suppose the lattice is uni- 
formly strained along the axis m. In a coordinate system 
aligned with this strain axis, we have 



A-Va 
A = | A- 1 / 2 
OA 



(6) 



where A is the distortion factor, i.e. the strained length 
normalized by the original length of the sample, which is 
related to the strain e by A = 1 + e. In a general coor- 
dinate system, the distortion tensor components become 



A Q/3 = \-^ 2 S, 



af3 ■ 



(a-a- 1 ^ 



y m/3. 



(7) 



The third tensor required for neoclassical rubber elas- 
ticity is the shape tensor £ at the time of crosslinking. 
For now, suppose the system is crosslinked in a totally 
disordered state, with no long-range or even local orien- 
tational order. In that case, a reasonable model for £q 1 
is the isotropic average of £ . This average gives the 
tensor components 



0,a/3 



"a/3, 



(8) 



(2£- 1 1 +£- 1 )/3. 



where a 

To determine the elastic term in the lattice Hamilto- 
nian, we substitute the tensor expressions ©, (Q), and 
© into the general formula of Eq. . In this substitu- 
tion, we note that det£ is constant because £ represents 
perfect local orientational order at a specific lattice site. 
Hence, the determinant term adds an unimportant con- 
stant to the Hamiltonian, and we can neglect it. After 
some algebra, the trace term leads to 



^clastic (fli , A) 



(9) 



(A 2 + 2A- 1 )- 7 (A 2 -A- 1 ) 



3 / \2 



In this expression, the first term is the classical elastic 
free energy for conventional isotropic elastomers, and 
the second term represents the anisotropy of liquid- 
crystalline elastomers. As expected, the second term 
shows a coupling between the elastic strain and the di- 
rector orientation. When the lattice is strained, each 
local director rii tends to align along the strain axis m, 
with an aligning potential that increases as the distor- 
tion A increases above 1. Conversely, when the directors 
are aligned, the lattice tends to extend along the average 
director. The strength of the coupling depends on the 
parameter 



2t]} - 2£7 1 

2£i x Vq l ' 



(10) 



which represents the difference in polymer step lengths 
parallel and perpendicular to the local director. This 
parameter expresses the anisotropy of the local mesogenic 
units, and controls how the local director interacts with 
the strain. It ranges from (in the isotropic case £\\ = £±) 
to 1 (in the maximally anisotropic limit £u 3> £±). 

We must now consider the possibility of symmetry- 
breaking fields acting on the elastomer. Symmetry- 
breaking fields can arise from two possible sources. The 
simplest possibility is a uniform stress a applied to the 
elastomer. Such a stress couples to the strain e, or equiv- 
alently to the distortion A = 1+e, and gives an additional 
contribution to the Hamiltonian of — crA for each lattice 
site. A more subtle possibility is a symmetry-breaking 
field quenched into the local shape tensor £o at the time 
of crosslinking. If the system has long-range order at the 
time of crosslinking, then £q is anisotropic with a sin- 
gle principal axis tiq at all lattice sites. If the system 
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has short-range order at the time of crosslinking, then £ l Q 
is anisotropic with a different principal axis n l at each 
lattice site i. 

In principle, we can incorporate the long- or short- 
range anisotropy of £ l Q into the model by writing a gen- 
eral expression for this tensor and substituting it into 
the trace formula of Eq. © . The detailed calculation is 
not algebraically tractable, but by symmetry we can see 



that the tensor components 



'0,al3 



must involve a com- 



bination of the isotropic tensor 6 a p and the anisotropic 
at site i. The anisotropic term acts as 



tensor nQ a n l p 
an effective field on the director rii, with a coupling of 
the form — (hi ■ rii) 2 . The direction of hi is the local 
quenched- in axis n l , and the magnitude of h t scales with 
the magnitude of the local quenched-in nematic order. 

By combining Eq. @, Eq. 10, and the effects of 
symmetry-breaking fields, we obtain the final expression 
for the lattice Hamiltonian 



F 



^ t Jij( n i 



E 



(id) 



M /\2 



(A 2 + 2A- X ) 



-a\ — (hi ■ rii 



(11) 



In this Hamiltonian, the statistical variables are the di- 
rector rii on each lattice site i and the overall lattice 
distortion A. We can assume that the distortion axis m 
is aligned with the principal axis of the director distribu- 
tion. 

This model of Eq. (|If \) can describe uniform elastomers 
or elastomers with quenched random-bond or random- 
field disorder. Uniform elastomers are represented by a 
bond strength Jy = J independent of position, and by 
a local field hi = 0. In this case, the system can have 
an isotropic-nematic transition with a transition temper- 
ature Tjjq that depends on J. Random-bond elastomers 
are represented by a bond strength JJy that depends on 
the position of the lattice sites i and j. We can regard this 
variation in the bond strength as a variation in the local 
Tin, which could be caused by chemical heterogeneity in 
the elastomer. Random-field elastomers are represented 
by a local field hi that varies randomly with position. 
This variation models randomness in the local orienta- 
tional order at the time of crosslinking, which gives het- 
erogeneous local stresses on the elastomer. We can now 
perform simulations to study the isotropic-nematic tran- 
sition in each scenario. These simulations are presented 
in the following section. 



III. SIMULATIONS 

We simulate the model of Eq. using the Monte 
Carlo method. We run the simulations on a 3D cubic 
lattice with periodic boundary conditions. We use a 3D 



rather than a 2D lattice, even though the simulations 
take longer in 3D, in order to avoid a 2D Kosterlitz- 
Thoulcss transition 20] . We normally use a lattice of size 
36 x 36 x 36. However, we have run a limited number of 
simulations on a larger lattice of size 48 x 48 x 48, for the 
uniform, random-bond, and random-field cases, and the 
results are generally consistent. In the simulations, we 
take the uniform or average value of the bond strength 
Jij to be I, and the shear modulus \i to be 1. These pa- 
rameters define a scale for the temperature. We let the 
anisotropy parameter 7 have its maximum value of 1, in 
order to see the greatest coupling between orientational 
order and elastic distortion. 

In each Monte Carlo step of the simulations, we at- 
tempt one local director rotation per lattice site and one 
change in the overall elastic distortion A. At each tem- 
perature, we equilibrate for 3000 Monte Carlo steps, and 
then collect data for 2000 Monte Carlo steps. This num- 
ber of steps is sufficient to reach equilibrium at all tem- 
peratures except for certain cases of hysteresis, which 
are discussed below. From the numerical data, we ex- 
tract two parameters as functions of temperature T: the 
elastic distortion A(T) and the orientational order pa- 
rameter S(T), which is defined as the largest eigenvalue 
of the tensor Q a p = (^ni^rii^ — ^5 a p), averaged over 
lattice sites i. That parameter represents the degree of 
ordering of the local directors along an average axis. 

We begin with the local directors in a disordered con- 
figuration, and then cycle the temperature downward and 
back upward, using the ending configuration at one tem- 
perature as the starting point for the next. This tem- 
perature cycling mimics the procedure in typical exper- 
iments, and provides an explicit test for hysteresis. For 
most parameter sets, we vary the temperature from 0.9 
to 0.8 and back in steps of 0.05, for a total of 41 runs 
over the temperature range. Each temperature cycle re- 
quires approximately 48 hours on a single processor of 
the Huinalu linux supercluster at the Maui High Perfor- 
mance Computing Center. 

Because the temperature cycle gives two runs at each 
temperature (except the lowest), we can assess whether 
each run shows a stable, metastable, or unstable state. 
To test for unstable states, we check whether the order 
parameter S has stabilized by fitting it as a linear func- 
tion of the Monte Carlo step number over the final 2000 
steps at each temperature. We identify a state as unsta- 
ble and remove it from our results if the absolute value of 
the slope exceeds a threshold. In practice, we find that 
the threshold of 3.2 x 10~ 5 eliminates at most two runs 
from each temperature cycle, one on cooling and one on 
heating. To test for stability vs. metastability, we com- 
pare the values of S for cooling and heating runs at the 
same temperature. If these values are within six stan- 
dard deviations of each other, we assume they represent 
the same stable state, so we average the two runs to ob- 
tain one data point with improved statistics. If not, we 
assume that one state is stable and the other metastable, 
so we report both in our results. 
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FIG. 2: (Color online) Plots of the orientational order param- 
eter S(T) and the elastic distortion A(T) for a homogeneous 
elastomer under several values of the applied uniform stress a. 
At zero stress, both plots show a first-order isotropic-nematic 
transition with a large hysteresis. As the stress increases, the 
first-order discontinuity decreases and then vanishes as the 
system passes through a mechanical critical point. Tempera- 
ture is in units of the bond strength J. 



Uniform Elastomers 



For an initial series of simulations, we consider uniform 
elastomers, with no randomness in the bonds (all Jjj = 1) 
and no random fields (all hi =0). The numerical results 
for these simulations are shown in Fig. 2. For zero ap- 
plied stress <t, the system has a first-order transition with 
hysteresis between the high-temperature isotropic phase 
and the low-temperature nematic phase. On cooling, the 
orientational order parameter S(T) jumps from 0.05 to 
0.62, and the lattice distortion A(T) jumps from 1.01 to 
1.22 (1% to 22% strain), at a scaled temperature of 0.81. 
On heating, S(T) jumps from 0.47 to 0.02, and A(T) 
jumps from 1.13 to 1.00 (13% to 0% strain), at a scaled 
temperature of 0.85. The large jumps in both of these 
order parameters, and the width of the hysteresis region, 
show the strong first-order character of the transition. 

When a symmetry-breaking stress is applied to the 
uniform elastomer, the phase transition changes drasti- 
cally. An applied stress increases both order parameters 
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FIG. 3: Schematic view of the block structure in the simu- 
lations with quenched disorder, (a) Random bond strength, 
(b) Random field orientations. 



S(T) and A(T) for all temperatures. As the stress be- 
comes larger, the transition temperature increases, the 
first-order jumps in the order parameters decrease, and 
the hysteresis region becomes narrower. At a critical 
value of the stress between 0.06 and 0.08, the first-order 
jumps vanish and the hysteresis goes away. Beyond that 
stress, the system shows a smooth supercritical evolution 
between the high-temperature disordered limit and the 
low-temperature ordered limit. As the stress continues 
to increase, this supercritical evolution becomes increas- 
ingly broad. This trend with increasing stress is consis- 
tent with the prediction of de Gennes based on symmetry 
considerations [TJ. It is analogous to the critical point 
in the liquid-gas transition under high pressure. 

The simulations show that the orientational order pa- 
rameter S(T) and the elastic strain e(T) = A(T) — 1 have 
roughly the same dependence on both temperature and 
applied stress. The linear scaling between S(T) and e(T) 
is consistent with the prediction based on symmetry con- 
siderations. 

We note that the smooth evolution of S(T) and e(T) 
beyond the mechanical critical point agrees with ex- 
periments on the isotropic-nematic transition in liquid- 
crystalline elastomers. However, as discussed in the 
Introduction, our previous paper found experimen- 
tal indications that the isotropic-nematic transition is 
smooth even if an elastomer is not under a supercrit- 
ical stress IT]]. For that reason, we need to look for 
other mechanisms to broaden this transition. Hence, we 
consider random-bond and random-field disorder in the 
following sections. 



B. Random Bonds 

To simulate disordered systems, we first consider elas- 
tomers with variations in the local bond strength. We 
let the parameters Jij of Eq. (|11|) be quenched random 
variables, which are fixed at the beginning of the sim- 
ulation and do not change in response to the statisti- 
cal evolution of the directors and the lattice distortion. 
We suppose that the variation of Jij occurs in blocks, as 
shown in Fig. 3(a). Within each block, Jy has a uniform 
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FIG. 4: (Color online) Orientational order parameter S(T) 
for simulations of elastomers with random-bond disorder, 
(a) Varying magnitude of disorder, with fixed block size 
12 x 12 x 12. (b) Varying block size of disorder, with fixed 
magnitude ±20%. As the disorder magnitude and block size 
increase, the transition is broadened into a smooth crossover 
between the isotropic and nematic limits. Temperature is in 
units of the average bond strength Jij. 



value from a Gaussian distribution with mean 1. From 
block to block, there are no correlations in . Both the 
width of the Gaussian distribution and the size of the 
blocks are parameters for the model, which are discussed 
below. Because the isotropic-nematic transition temper- 
ature depends on the bond strength, this model repre- 
sents an elastomer with blocks of different local transi- 
tion temperature Tjjv- This could occur if the elastomer 
is chemically heterogeneous, with different compositions 
in different local regions. 

In Fig. 4(a), we show the results for varying the mag- 
nitude of the disorder, i.e. the standard deviation of 
the Gaussian distribution of J^, for the fixed block size 
12 x 12 x 12. To save space, we show only the plots of 
the orientational order parameter S(T), because the cor- 
responding plots of the lattice distortion A(T) look quite 
similar. For weak disorder of ±5%, the system has a 
strong first-order isotropic-nematic transition, which is 
almost identical to the result for no disorder. Clearly the 
orientational order can average over this weak disorder 



strength. For a larger disorder of ±10%, there is still a 
first-order transition, but the first-order discontinuity in 
S(T) is smaller, and the width of the hysteresis region is 
greatly reduced. When the disorder reaches ±15%, there 
is no longer a first-order discontinuity nor a hysteresis 
loop. Instead, the material evolves rapidly but smoothly 
between the isotropic and nematic limits as a function of 
temperature. For an even larger disorder of ±20%, the 
transition becomes even broader, with a reduced slope in 
S(T). 

Figure 4(b) presents the results of varying the block 
size for a fixed disorder magnitude of ±20%. For the 
large block size 12 x 12 x 12 discussed above, or even for 
block size 6x6x6, the random-bond disorder causes 
a smooth evolution between the isotropic and nematic 
states. However, for a reduced block size of 3 x 3 x 3, 
there is a first-order isotropic-nematic transition, with 
a small first-order discontinuity and small hysteresis. If 
the block size is reduced to 1 x 1 x 1, i.e. each block is 
a single lattice site, then the system has a strong first- 
order transition. Although the disorder strength is quite 
large, the result for block size 1 x 1 x 1 is almost identical 
to the result for no disorder. Hence, reducing the block 
size is effectively equivalent to reducing the magnitude of 
disorder. The orientational order can average over small 
blocks of large disorder, just as it averages over large 
blocks of small disorder. 

The results of this section show that random-bond 
disorder can change the nature of the isotropic-nematic 
transition in liquid-crystalline elastomers, provided that 
the magnitude and length scale of the disorder are suf- 
ficiently large. If those conditions are satisfied, then the 
orientational order parameter S(T) undergoes a contin- 
uous change from a low value in the high-temperature 
isotropic limit to a large value in the low-temperature 
nematic limit. The lattice distortion A(T) goes through 
a corresponding smooth evolution. Thus, this type of 
disorder provides one mechanism to explain the experi- 
mental results. 



C. Random Fields 

As an alternative to random bonds, quenched disorder 
might affect liquid-crystalline elastomers through ran- 
dom fields coupling to the local directors. To simulate 
random-field effects, we let the fields hi of Eq. 1|11|) be 
quenched random variables, and let the bond strengths 
Jij be fixed at 1. We suppose that the fields hi have 
a fixed magnitude h and random orientation. As in the 
random-bond case, we suppose that the randomness oc- 
curs in blocks, as shown in Fig. 3(b). The random orien- 
tation is uniform at every site within a block, and it has 
no correlations from block to block. The magnitude h of 
the random field and the size of the blocks are thus two 
parameters for this model. Note that this model repre- 
sents an elastomer with different preferred orientations 
of the local director in different blocks. This could occur 
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FIG. 5: (Color online) Orientational order parameter S(T) 
for simulations of elastomers with random-field disorder, 
(a) Varying magnitude of random field, with fixed block size 
12 x 12 x 12. (b) Varying block size, with fixed random-field 
strength h = 0.4. As the random-field strength and block size 
increase, the transition is broadened and then the ordered ne- 
matic phase is destroyed. Temperature is in units of the bond 
strength J. 



if the crosslinking process quenches heterogeneous local 
stresses into the polymer network. 

Simulation results for several values of the random-field 
strength h at fixed block size 12 x 12 x 12 are presented 
in Fig. 5(a). For a small random field h = 0.2, the sys- 
tem has a first-order isotropic-nematic transition, which 
is fairly close to the result for no disorder. For a slightly 
larger random field h = 0.3, the magnitude of the first- 
order discontinuity and the width of the hysteresis region 
are both reduced. For h = 0.4, the first-order transition 
is much weaker, and the system is close to the smooth 
crossover between the isotropic and nematic phases seen 
in the previous two sections. However, for a larger ran- 
dom field h = 0.8, the behavior is quite different. Instead 
of broadening the isotropic-nematic transition, the ran- 
dom field simply destroys the long-range nematic order 
in an athermal way. In that high-field limit, each local 
director is just aligned with its local random field, giving 
a slight residual order parameter S(T) that is approxi- 
mately independent of temperature. 



Figure 5(b) shows the simulation results for several val- 
ues of the block size at the fixed random-field strength 
h = 0.4. As in the random-bond case, varying the 
block size has the same effect as varying the random- 
field strength. For small block size, the system has a 
first-order transition that is very close to the result for no 
disorder. This behavior shows that the orientational or- 
der averages over small blocks of strong random field. For 
larger block size, the phase transition gradually changes 
toward a smooth crossover between the isotropic and ne- 
matic limits. Analogous simulations for h = 0.8 (not 
shown) demonstrate that increasing the block size takes 
the system from a strong first-order transition toward 
a limit in which orientational order is destroyed for all 
temperature. 

Overall, the simulations presented in this section 
show that random fields can have two distinct effects 
on the isotropic-nematic transition in liquid-crystalline 
elastomers. Low random-field disorder broadens the 
isotropic-nematic transition, but high random-field dis- 
order destroys the nematic phase in a temperature- 
independent way. The results for low random-field disor- 
der are consistent with experiments on liquid-crystalline 
elastomers, but the results for high random-field disorder 
show a very different type of behavior. 



IV. DISCUSSION 

In this paper, we have investigated three possibile 
mechanisms to explain the broadening of the isotropic- 
nematic transition in liquid-crystalline elastomers. The 
first mechanism is a stress that couples to the strain and 
hence to the nematic order. Our simulations show that 
a small stress reduces the first-order discontinuity in the 
isotropic-nematic transition, and a critical stress causes 
this discontinuity to vanish. Beyond that point, the elas- 
tomer has a smooth crossover between the isotropic and 
nematic limits as a function of temperature. Hence, a 
supercritical stress gives the same general trend as the 
experiments. However, our previous paper found evi- 
dence that the smooth isotropic-nematic transition does 
not require a supercritical stress ^lj. I n particular, a 
smooth transition is seen even in elastomers crosslinked 
in the isotropic phase, even under the minimum applied 
stress required for the experiment, conditions in which a 
supercritical stress is unlikely to occur. 

An alternative possibility to explain the broadened 
transition is quenched disorder in the elastomer. One 
specific mechanism to generate quenched disorder is 
chemical heterogeneity, which can be represented by ran- 
dom bonds in a lattice model. Our simulations show that 
random-bond disorder can broaden the isotropic-nematic 
transition into a smooth crossover, if the magnitude and 
length scale of the disorder are large enough. This 
numerical result is consistent with theoretical work on 
generic random-bond systems by Imry and Wortis |2l| . 
which argued that weak random-bond disorder should 
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reduce the first-order discontinuity in a transition, and 
larger disorder should eliminate the discontinuity com- 
pletely. It is also consistent with other simulation stud- 
ies of the isotropic-nematic transition in systems of small 
molecules with quenched random impurities |'2'2t |23| . 
Thus, random-bond disorder provides a plausible mecha- 
nism to explain experimental results on liquid-crystalline 
elastomers. 

Another type of quenched disorder is heterogeneous lo- 
cal stresses, which can be represented by random fields in 
a lattice model. In general, random fields have stronger 
effects on ordered phases and phase transitions than ran- 
dom bonds. In our simulations, we find that low random- 
field disorder broadens the isotropic-nematic transition, 
consistent with the experiments, but high random-field 
disorder destroys the nematic phase over all tempera- 
tures. 

Our random-field simulation results are surprising in 
comparison with theoretical expectations for random- 
field systems. In classic work on quenched disorder, Imry 
and Ma [24| showed that arbitrarily small random fields 
should destroy long-range order in a discrete order pa- 
rameter (such as an Ising model) for spatial dimension 
less than 2, and destroy long-range order in a continuous 
order parameter for spatial dimension less than 4. In our 
case, the nematic order parameter is continuous, and the 
spatial dimension is 3. Hence, the Imry-Ma argument 
implies that random fields should destroy nematic order. 
The simulations of high random fields show this effect, 
but the simulations of lower random fields show a broad 
isotropic-nematic transition, which is a very different ef- 
fect. 

One might ask whether this behavior results from an 
Imry-Ma domain size that is large compared with the 
size of the simulation cell. To answer that question, we 
can estimate the Imry-Ma domain size. In the Imry-Ma 
argument, domains form with a characteristic size £ such 
that the boundary energy equals the field energy. In our 
3D system of continuous directors, the boundary energy 
is of order J£. To estimate the field energy, recall that 
we are simulating blocks of sites with the same random 
field, as shown in Fig. 3. Let b be the linear size of a 
block. The field energy for a single block is of order 
Sbiock (hi ' n i) 2 ~ h 2 b 3 , and the number of blocks per 
domain is of order (£/6) 3 . As a result, the field energy for 
a domain is of order (/i 2 6 3 )(£/6) 3 / 2 w h 2 b 3 / 2 £ 3 / 2 . This 
argument implies that increasing the block size causes the 
random field to be more effective, as is seen in the sim- 
ulations. Comparing the boundary energy with the field 
energy gives an Imry-Ma domain size of £ « J 2 /(h 4 b 3 ). 
For the simulated values J = 1, h = 0.4, and b = 12, this 
size is less than I lattice unit, much less than the system 
size. Thus, the behavior in our simulations does not arise 
from a large Imry-Ma domain size. 

We suggest that the new behavior in our simulations 
arises from the coupling between the local directors and 
global elastic strain variable. The Imry-Ma prediction is 
based on an analysis of the energetics of local ordered 



domains. This analysis assumes a short-range interac- 
tion in the order parameter. By contrast, in our model 
for liquid-crystalline elastomers, the local director at any 
lattice site interacts with the global elastic strain, which 
in turn interacts with the local director at every other lat- 
tice site. Hence, the elastic strain mediates an effective 
long-range interaction between local directors on differ- 
ent sites. This changes the assumptions in the Imry-Ma 
theory, and hence allows a broad isotropic-nematic tran- 
sition over a range of random-field strength. 

For a specific numerical test of this suggestion, we per- 
form simulations of a simplified model without the globcl 
elastic strain variable. The Hamiltonian for this model 
consists only of the interaction of Eq. plus the random 
fields acting on the local directors. In this case, the in- 
teraction is purely short-range, so the Imry-Ma argument 
should apply. Indeed, the simulations show that random 
fields simply destroy the nematic order and do not in- 
duce a broad isotropic-nematic transition. This confirms 
the concept that a broad transition is a new effect arising 
from a strain-mediated long-range interaction. 

In a realistic system, the elastic strain-mediated in- 
teraction is not infinite-T&nge, as in our model. How- 
ever, elastic interactions do have a power-law form, and 
hence they can have long-range effects. In a recent 
renormalization-group study, Xing and Radzihovsky [l^ 
have assessed the effects of elastic interactions on liquid- 
crystalline elastomers with random fields. They find that 
elastic interactions cause the nematic order to be robust 
against the disordering effect of random fields. Through a 
power-law expansion about spatial dimensionality 5, they 
estimate that nematic order can be stable down to a crit- 
ical dimension well below 3. This is apparently the same 
stabilization that we see numerically. Thus, our simula- 
tion shows the consequence of this stabilization for the 
temperature-dependent isotropic-nematic transition. It 
is interesting to note that similar considerations of disor- 
der and long-range elasticity have been seen in models for 
magnetic pha se transitions in colossal magnetoresistancc 
materials |2f|. 

In conclusion, we have developed a lattice model of 
the isotropic-nematic transition in liquid-crystalline elas- 
tomers. The model considers a local directors cou- 
pled to a global elastic distortion variable, and allows 
both random-bond and random-field disorder. Through 
Monte Carlo simulations of this model, we find that a 
uniform elastomer has a mechanical critical point, that 
both random-bond disorder and low random-field disor- 
der broaden the isotropic-nematic transition, and that 
high random-field disorder destroys the nematic phase. 
The model therefore confirms that the width of the 
isotropic-nematic transition can be controlled by hetero- 
geneity in liquid-crystalline elastomers. 
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